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ABSTRACT 

Context. During solar flares a large number of charged particles are accelerated to high energies, but the exact mechanism responsible 
for this is, so far, still unclear. Acceleration in collapsing magnetic traps is one of the mechanisms proposed. 

Aims. In the present paper we want to extend previous 2D models for collapsing magnetic traps to 3D models and to 2D models with 
shear flow. 

Methods. We use analytic solutions of the kinematic magnetohydrodynamic (MHD) equations to construct the models. Particle orbits 
are calculated using the guiding centre approximation. 

Results. We present a general theoretical framework for constructing kinematic MHD models of collapsing magnetic traps in 3D and 
in 2D with shear flow. A few illustrative examples of collapsing trap models are presented, together with some preliminary studies 
of particle orbits. For these example orbits, the energy increases roughly by a factor of 5 or 6, which is consistent with the energy 
increase found in previous 2D models. 

Key words. Sun: corona - Sun: flares - Sun: activity - Sun: magnetic fields - Sun: X-rays, gamma rays 



ON 

o 
o 

(N 
+-> 

o 
O 
o 

<N 

5? 
43 

6 

c3 1. Introduction 



> 

in 
oo 
m 

O 

On 
O 



One of the main features of solar flares is the acceleration to 
high energies of a substantial number of charged particles within 
a short period of time. The explanation of how this happens is 
one of the most important open questions in solar physics. There 
is general agreement that the energy released in solar flares is 
previously stored in the magnetic field, but the exact physical 
mechanisms by which this energy is released and converted into 
bulk flow energy, thermal energy, non-thermal en ergy and ra- 
diation energy are still a matter of discussion (e.g. [M iller et al. 
l~997t lAschwanderJ|2002t lNeukirchll200H iNeukirch et alj |2007t 



magnetic trap (CMT from now on) the kinetic energy of the par- 
ticles could increase due to the betatron effect, as the magnetic 
field strength in the CMT increases, and due to first-order Fermi 
acceleration, as the distance between the mirror points of particle 
orbits decreases due to the shortening of the field lines. There is 
also some observational evidence of post-flare field lines relax- 
ation (field line shrink age) from Yohkoh ( e.g. iForbes & Actonl 
1996) and Hinode (e.g. lReeves et alj|20 08b) observations. 

Various fundamental properties of the particle acceler- 
ation process in C MTs have been investigated by Somov 



Krucker et alj|200a lAschwand en 2009). Using observations of 



non-thermal high-energy (hard X-ray and y-ray) radiation, it is 
estimated that a large fraction of the released magnetic energy 
(up to the order of 50 %) is converted into non-thermal energy i n 
the form of high energy particles (e.g. Ems lie et alj2004 12005). 

A variety of possible particle acceleration mechanisms have 
been suggested including direct acceleration in the parallel elec- 
tric field associated with the reconnection process, stochas- 
tic acceleration by turbulence and/or wave-particle resonance, 
shock acceleration or acceleration in the inducti ve electric 
field of the reconfiguring magnetic field (see e .g. Miller etaT] 
19971: lAschwanden! 20021 : lNeukirchll200H INeukirch et alj|2007t 
Krucker et alj 12008, for a detailed discussion and further ref- 



erences). So far, none of the proposed mechanisms can ex- 
plain the high-energy particle fluxes within the framework of 
the standard solar flare thick target model. This has recently 
prompted suggestions of alternative accel eration scenarios (e.g. 
iFletcher & Hudson 2008; iBim et afll2009l). 



Somovl d 19921) and ISomov & Kosugil d 19971) suggested that 



the reconfiguration of the magnetic field during a flare could 
contribute to the acceleration of particles. Due to the geometry 
of the magnetic field charged particles could be trapped while 
the magnetic field lines relax dynamically. In such a collapsing 



and co-workers (e.g. iBogachev & Somovl 1200 U 120 05, 2009; 
iKovalev & Somovl |2002t l2QQ3allbh ISomov & Bogachevl 120031) . 
including the relative efficiencies of betatron and Fermi accel- 
eration, the effect of collisions, the role of velocity anisotropies 
and the evolution of the energy distribution function in a 
CMT. In all cases a basi c model for CMTs has been used. 
Karli ckv & Kosugil d2004l) also investigated particle accelera- 
tion, plasma heating and the resulting X-ray emission using a 
simple C MT model and a s i mplifi ed equation of motion for the 
particles. Karlickv & Barta (2006) used CMT-like electromag- 
netic fields taken from an MHD simulation of a reconnecting 
current sheet to investigate acceleration using test particle cal- 
culations with a view to explain hard X-ray loop-top sources. 
Some indication that CMTs might be relevant for X -ray loop top 
sources has been provided by Veronig et al. (2006)) . A very sim- 
ple tim e-dependent trap model was also used by lAschwanden! 

(2004) to explain the pulsed time profile of energetic particle in- 
jection during flares. 

A general theoretical framework for more detailed analytical 
CMT models based on kinematic MHD, i.e. with given bulk flow 

I I 

profile, in Cartesian coordinates was presented by Giuliani et al. 

(2005) for 2D and 2.5D magnetic fields, but excluding flow in 
the invariant direction. Some examples of model CMTs were 
given together with a calculation of a particle orbit based on 
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non-relativistic guiding centre theory (see e.g. iNorthropll 1963b . 
It was found that in the models studied the curvature drift and 
the gradient-B drift play an important role in the acceleration 
process. Similar findings have also been made, albeit in systems 
of a much smaller length, in the investigation of particle accel- 
eration in particle-i n-cell simulations of collisionless magnetic 
reconnection (e.g. Hoshino et al. 2001b. 

The advantage of kinematic MHD models compared to e.g. 
MHD simulations is that they allow us to obtain analytical ex- 
pressions for the electromagnetic fields of the CMT This makes 
the integration of particle orbits more accurate, because there 
is no need for interpolation of the fields between grid points. 
Furthermore the investigation of different model features is pos- 
sible in an easy way by varying model parameters. The major 
disadvantage of kinematic MHD models is their lack of self- 
consistency, but this is not too critical for the purpose of test 
particle calculations. 

Particle acceleration through rapid reconfiguration of the 
magnetic field has also been identified as one of the mechanisms 
for particle energization duri ng magnetospheric substorms (e.g. 
iBirn et alj[l997l Il998l |2004|) . During a substorm the stretched 
magnetic field of the magnetotail reconnects, leading to a so- 
called dipolarisation of the near-Earth tail, which is in principle 
very similar to the evolution of the magnetic field in a CMT as- 
sociated with a solar flare. A general comparison of flare and 
substorm/magnetotail ph enomena based on observations has re- 
cently been presented bv lReeves et alj (2008a). 

The purpose of the present paper is to extend the theo- 
retical framework fo r kinematic MHD CMT models given by 
iGiuliani et al.l d2005l) to 2.5D models with flow in the invariant 
direction and to fully three-dimensional models. This is neces- 
sary for a number of reasons: 



calculation of Euler potentials and uses the magnetic field 
directly. 



The paper is organised as follows. In Sect. 12.11 we briefly 
summarise the present state of the kinematic MHD theory of 
CMTs, before presenting its extensions to 2.5D with shear flow 
and to 3D. In Sect. [3]a couple of illustrative examples of CMT 
models based on the new theoretical descriptions are shown, fol- 
lowed by examples of test particle calculations in Sect. [4] We 
conclude the paper in Sect. [5] with a summary and conclusions. 
Appendix A gives more detail of the calculation of the 3D field 
using Euler Potentials. 



2. Basic Theory 

The CMT is assumed to form outside the nonideal reconnection 
region, so the ideal kinematic MHD equations may be used to 
describe the evolution of the electromagnetic field 



E + vxB = 0, 

i)B „ „ 

— = -VxE, 
dt 

V B = 0, 



(1) 
(2) 
(3) 



with the MHD velocity v assumed to be given as a function of 
space and time. We will also make occasional use of the ideal 
induction equation 



Y-Vx(vxfi), 
at 

which results from combining Eqs. (fl} and (fJJ. 



(4) 



1. The theory of kinem atic MHD CMTs as developed so far by 
Giuli arfiet al.l (2005) only allows for a magnetic field com- 
ponent in the invariant direction, but not for a component of 
the flow velocity in this direction. Without this component of 
the flow velocity the magnetic field component in the invari- 
ant direction can only increase in a CMT due to magnetic 
flux conservation. It is, however, to be expected that during 
a flare magnetic shear will be reduced rather than increased 
and therefore the introduction of a component of the flow ve- 
locity in the invariant direction is a necessary extension to be 
able to make the 2.5D models more realistic. 

2. In the 2D cases investigated by IGiuliani et al.l d2005l) the ac- 
celeration due to curvature and gradient-B drift occurs in the 
invariant direction. This is due to the fact that the particles 
gain energy while moving parallel or anti-parallel (in the 
case of electrons) to the inductive electric field which in a 
2D trap is in the invariant direction. Due to the spatial sym- 
metry the electric field does not vary in this direction and 
this will have an influence on the acceleration process. It is 
therefore important to investigate the differences of the ac- 
celeration process between 2D models and non-symmetric 
3D CMT models in t he future. 

3. IGiuliani et al.l (2005) have already discussed a possible way 
of extending the 2D theory to three dimensions using Euler 
potentials. While Euler potentials allow a relatively straight- 
forward extension of the theory to 3D by simple analogy to 
the 2D case, they are not easy to use in the modelling pro- 
cess, which is already intrinsically more difficult in three di- 
mensions. We therefore present in this paper an extension 
to the theory which makes it possible to avoid the explicit 



2.1. Kinematic MHD Models of CMTs in 2.5D without shear 
flow 

We start by giving a brief overview of the translationally in- 
variant 2.5D kinema tic MHD theory of CMTs developed by 
IGiuliani et alj (I2005I) . This does not include a velocity compo- 
nent in the invariant direction. In t he following we will use the 
same coordinate system as used bv lGiuliani et alj d2005l) . i.e. all 
physical quantities depend only upon x and y, with x being the 
coordinate parallel to the solar surface (photosphere) and y being 
the height above the solar surface. The invariant direction is the 
z-direction. 

For the cases with spatial symmetry it is useful to write the 
magnetic field as 

B = B p + B z e z = VA x e z + B z e z , (5) 

where A(x, y, t) is the flux function, B p = 
(B x (x,y, t),B y (x,y,t),0) and B,(x,y,t) the z-com ponent of the 
magne tic field. An important assumption made by Giulian i~et al.l 
d2005l) is that there should be no flow in the invariant direction, 
i.e. 



v 2 (x,y,t) = (v x (x,y,t),v y (x,y,t),0). 



(6) 



As we will make use of this particular velocity field later on, we 
use the index 2 to distinguish it from the full velocity field with 
non-zero v z . Using an appropriate gauge for A, the z component 
of Ohm's law ([TJ gives 



dA 8A _ t n 
— = — + v 2 ■ VA = 

dt dt 



(7) 
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for the time evolution of the flux function A. For the time evo- 
lution of B z it is better to use the z-component of the induction 
equation (|4), 

^ + V ■ (v 2 B z ) = 0. (8) 
ot 

Equations © and dHJ simply express the conservation of mag- 
netic flux. In the case with vanishing shear velocity (v, = 0) 
the magnetic flux J B z dxdy is conserved ind ependently. To solve 
Eqs. and © for A(x,y, f) and B z (x,y, f), lGiuliani et alj d2005l) 
prescribe a time-dependent transformation between Lagrangian 
coordinates X, Y and Eulerian coordinates x, y: 

X = X(x,y,t), Y=Y(x,y,t), (9) 

instead of a time-dependent velocity field v x (x,y, t), v y (x,y, t). 
The velocity field can be determined easily from th e trans forma- 
tion equations (see Eqs. (23)-(26) of Giulia nTet al.l d2005l) ). 

The solution for the magnetic flux function A(x, y, t) is then 
trivially given by 

A(x,y,t) = A (X(x,y,t),Y(x,y,f)\ (10) 

where Aq(X, Y) is the flux function at some reference time t = to. 
The B x - and B v -components of the magnetic field can be calcu- 
lated from Eq. <(5j by differentiation. 

Equation dHJ has the form of a continuity equation for B z 
with the solution 

lv ,„ .. a .. an 

fo and \J\ 



2.3. Extension to 3D 



B z {x,y,t) = r'B 0z (X(x,y,t),Y(x,y,t)), 

where Bo z {X, Y) is again the B z at a reference time t 
is the Jacobian determinant of the transformation between the 
Lagrangian and Eulerian coordinates, here written as 



(12) 



dX dY dY dX 
dx dy dx dy 

The Jacobian determinant basically expresses the deformation of 
infinitesimal area elements , i.e. a change of cross section, in the 
jc-y-plane during the time evolution of the system. Because the 
magnetic flux associated with B z is conserved independently in 
the case discussed in this section, any decrease in area must be 
compensated by a matching increase in B z and vice versa. 

Finally, the electric field can be determined from Ohm's law 
(Q]l once the velocity field v and the magnetic field B are known. 



2.2. Extension to 2.5D with shear flow 

To allow the effect of shearing and also de-shearing of the mag- 
netic field to be taken into account it is necessary to have a non- 
zero v z (x, y, f). The basic effect of a non-zero v z is to add a source 
term to equation © 



dB 



+ V ■ (v 2 B z ) = V ■ (v z B p ) . 



(13) 



The source term on the right-hand-side of Eq. (fT3l destroys the 
separate conservation of magnetic flux in the z-direction, be- 
cause a non-zero v z allows B x and B y to by turned into B z and 
vice versa. In addition to the transformation equations for the x- 
and y-coordinates one has to add a transformation equation for 
the z-coordinate of the form 

Z = z + Z(x,y,t). (14) 

The general solution for the flux function remains the same, but 
the solution for B z becomes more complicated. As it is much 
easier to deduce the solution for B z as a special case from the 3D 
case discussed next, we will give the appropriate expressions for 
B z and the velocity field after discussing the general theory for 
three dimensions. 



As already pointed out by Giulian i~et alj d2005l) . one can in prin- 
ciple use a similar approach as for 2D to generalise the theory to 
3D. Instead of writing the magnetic field in terms of a flux func- 
tion A we u se Euler Potentials to satisfy the solenoidal condition 
© (see e.g. lSte^fl970l[T987h : 



B = VaX VB. 



(15) 



When using Euler potentials one has to assume that the mag- 
netic topology of the CMT is sufficiently simple to allow the 
global existence of a set of Euler p otentials satis fying Eq. (15[ 
for all positions and times (see e.g. iMoffatdfl 978l for a discus- 
sion). Using Euler potenti als in an app ropriate gauge, Ohm's law 
(fl]l can be written as (e.g. Sternl ll97 



ppr 
0) 



da _ 

— + v Va = 0, 

dt 



and the solutions of Eqs. ( TT6b and dTTb are given by 

a(x,t) - a(X(x,t)), 
B(x,t) = j8(X(x,t)). 



(16) 
(17) 



(18) 
(19) 



where, as in the 2D solution a (X) and ]3(X) are the Euler po- 
tentials at a reference time t = foQ As in the 2D case a trans- 
formation between Eulerian (jc) coordinates and Lagrangian (X) 
coordinates is assumed as given in the form 



X = X(x,y,z,t), 



(20) 



where we have combined the transformation equations for the 
three coordinates (X = X(x, y, z, t), Y = Y(x, y, z, t), Z - 
Z(x, y, z, 0) into a vector X = (X, Y, Z) for ease of reference. For 
completeness, the full derivation of the expression for the mag- 
netic field using Eqs. (fT8l and ( fT9l is shown in Appendix|A] The 
result is given by the equations 



By = 



dX dX 

dz dx 

dX dX 

dx dy 



B (X), 



B z = \-r-x—\-B Q (X) 



(21) 
(22) 
(23) 



It is important to note that this result is expressed completely 
in terms of derivatives of the transformation equations and the 
magnetic field at the reference time t — to 



B = Bo (x,y,z), 



(24) 



i.e. no reference to Euler potentials has to be made when mod- 
elling CMTs in 3D. This is no surprise as the same result 
can also be fo und without the use of Euler potentials (see e.g. 
Moffatn [l978L p. 44), but using Euler potentials makes the tran- 
sition from 2D to 3D a bit more obvious. While Euler poten- 
tials are often very useful for gaining better theoretical insigh t 



tials are orten very userul tor gaining better theoretical insight 
(e.g. IStern|[T970l iHesse & Schindlerill988t iHesse et al.ll2005h . 
they are usually quite d i fficult to use for modelling purpos es (e.g. 
iPlatt & Neukirch! 1994tlRomeou & Neukirchll999ll2002l) . Also, 
due to this result the conditions for the global existence of Euler 



1 For example, iGiuliani e t al. ( 2005) use the final time as reference 
time. 
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potentials do not apply for the modelling of 3D CMTs and the 
modelling process is thus much less restrictive. It is therefore 
very beneficial to have a formulation which is based purely on 
the magnetic field at the reference time and on the transforma- 
tion equation (l20t . both of which we are free to choose. 

From the transformation equation d20l ). one can calculate the 
flow velocity by using that 



_ = _ +(v . V )X = 0, 



(25) 



from which one can calculate the velocity v by inversion of the 
non-singular 3x3 matrix VX, giving 



v(x,y,z,t) = -(VX)- 



ax 

dt ' 



(26) 



We refrain from giving the complete explicit form of the velocity 
field here, as it is rather lengthy and not too instructive. Finally, 
knowledge of the flow velocity and the magnetic field allows the 
calculation of the electric field from Ohm's law (Q~|). 

2.4. Derivation of the 2.5D case with shear flow formulae 
from the 3D case 

We will now come back to the 2.5D case with shear flow. The 
transformation equation ( TBI for the z-coordinate implies that 



8X 
~dz 



= (0,0,1). 



(27) 



Using Eq. d23l , the z-component of the magnetic field for the 
2.5D case with shear is given by 



B z (x,y,t) = 



dY dZ dZ dY 
dx dy dx dy 
dZ dX dX dZ 
dx dy dx dy 



BoAX) 



B y(X) 



dXdY_dYdX\ BQ(X) 
dx dy dx dy J 



(28) 



The last term of Eq. ((28) is identical to the 2.5D solution for 
B z without shear flow given in Eq. ( fTTT i. The other two terms 
represent the extra possibility of turning B x or B v flux into B : 
flux and vice versa. 

The velocity field can be determined by using the transfor- 
mation equations for the 2.5D case with shear flow in Eq. ( f25l >. 
This gives the components of the velocity as 



_ / dXdY dXdY\ldXdY dY dX 

x \ dt dy dy dt J\ dx dy dx dy 

I dXdY + dXdY\ldXdY dY dX\ , 

1 \ dx dt dt dx)\ dx dy dx dy , 



(29) 
(30) 



v, = - 



<9Z _ 
~dt 

8Zl 
dy\ 



32 

dx 

dXdY 
dx dt 



dX dY dX dY\ 
~dt"dy + ~dy~~dt) 



dXdY 
'dtlh 



dXdY 
dx dy 



dY dX 

dx dy 



-l 



3. Illustrative Examples of Collapsing Trap Models 

In the following we shall discuss some simple illustrative exam- 
ples of CMTs in 2.5D with shear flow and in 3D. Our main pur- 
pose here is to compare some of t he features of t hese e xtended 
models with the results found by iGiuliani et alj d2005l) for 2D 
m odels. Therefore, we s hall use one of the transformations used 
bv lGiuliani et alj (1200 5). We do not suggest that these examples 
can be regarded as realistic models of a flare, but they offer some 
insight into the basic features of 2D and 3D collapsing trap mod- 
els. 



3.1. An illustrative example for a 2.5D CMT model with shear 
flow 

We first add a shea r flow to the main example presented in 
IGiuliani et alj (I2005I) . Therefore, as in their paper, the 2D mag- 
netic field is generated using the flux function 



y + d/L\ 

A Q = cj arctan — - 

\x + 1/2/ 



jyo + d/L 

c\ arctan 

U-l/2 



(32) 



which represents a loop between two line currents at xq - +L/2, 
i.e. separated by a distance L and placed at a distance yo = —d 
below the photosphere. The magnetic field generated by the flux 
function ( 1321 is potential if regarded as a function of xq and yo- 
This potential field is the final field to which the CMT relaxes as 
t — » oo. In the model presented in this section, the magnetic field 
in the z-direction is set to zero, B, = 0, as t — > oo. 

At other times, the magnetic field will be non-potential and 
we will choose a coordinate transformation which gives an ini- 
tially sheared magnetic field, i.e. with B 7 + . To ensure conti- 
nuity from the model of IGiuliani et alj (120051) to our model the 
transformations of the x- and y-coordinates are the same as in 
their paper, i.e. 



yo = (at) 6 In 



1 + 



(atf 



1 + tanh[(y - L v IL)a\\ 



1 + tanh[(y - L v /L)ai] 
+ < : > y- 



(33) 



(34) 



This transformation basically stretches the magnetic field in the 
y-direction above a height given by L v /L, where the transition 
between unstretched and stretched field is contro lled by the pa- 
ramete r fli. We use the same parameter values as IGiuliani et al.l 
(120051) . namely a = 0.4, b = 1.0, L,/L = 1 and a\ = 0.9. For 
simplicity, the transformation depends on time only through the 
function yo(y, t). This time-dependence lets the field collapse to 
the final field described above as for t — > oo, y Q tends to y. Other 
important features of the transformation are that the foot points 
of magnetic field lines do not move during the collapse as for 
y — we have yo = for all t . 

Th e important difference to the model used bv lGiuliani et al.l 
(120051) is that we introduce an additional transformation for the 
z-coordinate giving rise to a shearing flow as discussed above. 
The transformation for the z-coordinate is chosen as 



— — - — — — — - — I . (31) zo =z + 6\yo(y,t)-y] 



°2.5D + X 



(35) 



Again, the electric field can be calculated from Ohm's law ([]]), 
once the velocity field and the magnetic field are known, but due 
to the complexity of the expressions we do not state them here 
explicitly. 



where 5 and A2.5D are parameters that are explained later. This 
transformation induces an x- and y-dependent shear motion. The 
reasoning behind choosing the transformation as given is as fol- 
lows: 
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Fig. 1. Field lines in the example 2D case with shear flow. 
Lengths are normalised to L — IQMm. The left plot shows the 
magnetic field at 1.05s, the right plot shows it at 50.8s. The col- 
lapse of the field lines in the y-direction is obvious. Note the 
difference in scale between the x-z-plane and the y-direction. 




Fig. 2. Top views of the field shown in Fig.Q] again at t= 1.05s 
and 50.8s. These plot show more clearly how the magnetic field 
unshears. 



1. the shear flow should be anti-symmetric with respect to x 
and vanish as |x| — > oo, which is achieved in a simple way 
by the x-dependence of the transformation and controlled by 
the parameter 02.50; 

2. the shear flow should vanish at the photosphere (no foot point 
motion) and be of noticeable strength only in the stretched 
area of the magnetic field and it should also vanish as t — > oo; 
this is achieved in a simple way by the y-dependence of the 
transformation; 

3. we should be able to control the magnitude of the shear flow, 
which is done by the parameter 6. 

An example of the effect of the transformation on the mag- 
netic field is shown in Figs.[T]to[3]for two different times (in our 
normalisation these are t — 1.05s and t = 50.8s. In these plots 
the values 6-1 and «2.5D = 1 have been used. The initial shear 
and stretching as well as the collapse and unshearing of the field 
are obvious when comparing the plots of the magnetic field for 
the two different times. 

3.2. An illustrative example for a 3D CMT model 

To generate an example model for a 3D CMT, we use two 
magnetic point sources placed at positions (-L/2, -d, 0) and 
(L/2, -d, 0), so the sources are located underneath the x-axis at 
depth -d under the photosphere (y = 0), and they are separated 
by a distance L. The potential magnetic field generated by these 




Fig.3. Side views of the field shown in Fig. Q]at t=1.05s and 
50.8s. One can clearly see the collapse of the field lines in the 
CMT. 



sources is then given by 

[(x + §) e.v + (y + d) e y + zo e z ] 

Bo = c\ 

[(x + L/2) 2 + (yo + d) 2 + zg] 

\(xo - 3 ) e, + (yo + d) e y + z e z ] 

_Cl 7 ^72? (36) 

[(xo - L/2) 2 + (y + d) 2 + zlY 

The value of c\ is chosen so that the maximum value of the mag- 
netic field on the photosphere is around 0.01 T (100 G). We 
choose c\ to be negative so that the magnetic polarity is nega- 
tive for xo positive. As in the previous 2D case, we use d — L 
and our standard normalisation L — 10 7 m. As in the 2D case 
with shear flow, this potential field is the final field to which the 
CMT relaxes as t — * oo. It can be consi dered as a 3D generalisa- 
tion of the 2D magnetic field used by iGiuliani et al.l (120051) and 
in the present paper in Sect. 13.11 

For this 3D example we choose a transformation which ini- 
tially twists the field lines around the y-axis above a given height 
and for a given distance fro m the y-axis, as well a s stretching 
them in the y-direction as in IGiuliani et al.l d2005l) . The time- 
dependence of the transformation then untwists the field while 
it relaxes. To achieve this feature we now transform the x- 
coordinate as well as the z-coordinate, while keeping the trans- 
formation for y as given in Eq. ( [34-b to make t his illustrative 
examp le more easily comparable to the work by Giulianj^Tal] 
(2005) and the 2.5D case with shear flow described above. 

The general structure of the x- and z-transformations is simi- 
lar to the 2.5D case, with the difference that the x-transformation 
now also depends on z, while the z-transformation depends on x 
as follows: 

z 

xo = x-8(y (y,i)-y)- i -, (37) 

a^ D + x L + z l 

Zo = z + S(y ( y ,t)-y)- — -. (38) 

+ x z + z l 

The y-dependence has the same effect as for the 2.5D case with 
shear flow, whereas parameters d and 030 control the amount of 
twist and the distance from the y-axis for which there is twisting. 
The form of the transformation ensures that field lines which 
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Fig. 4. Field lines for the 3D example case. Lengths are nor- 
malised to L — IQMm. Left image shows the trap at 1.05s, right 
shows once it has collapsed at 50.8s. We point out that there is 
difference in scale between the y-direction, extending from to 
10 L in the plot and the x-z-plane which extends between -2 L 
and 2 L in both directions. 




Fig. 5. Top views of the field shown in Fig. [4] again at t= 1.05s 
and 50.8s. These plot show more clearly how the magnetic field 
unshears. 



pass through the region where the transformation deviates no- 
ticeably from the identity transformation are twisted in the coun- 
terclockwise direction apart from being stretched in the y direc- 
tion. 

We show an example with parameter values of 6 — 0.001 
and «3d = 1 in Figs. |4j [5] and [6] Figure [4] shows the how field 
lines relax between the initial time (1.05 s in the normalisation 
used for this example) and a later time (50.8 s). Apart from the 
collapse built into the example by the ^-transformation (see Fig. 
[6]) we can clearly see the effect of field line twisting through the 
x- and z-transformations, in particular in Fig. [5] 

4. Test Particle Orbits 

We present a couple of example calculations of particle orbits 
for the CMT models presented abo ve to compare them to the 
case studied in lGiuliani et alj d2005l) . As the gyroperiod and gy- 
roradius of electrons are far smaller than the typical time and 
length scales of the collapsing traps we can use guiding centre 
theory to determine the particle trajectories (e.g. Northrop 19631; 
IGiuliani et al]|2005l) . 

Initial conditions for the test particles for both the 2D exam- 
ple with shear flow an d 3D example were ch osen to be compara- 
ble to those studied in Giuli anTet al.l (H005), i.e. we have set the 
particles to start at the point x = 0.1, y = 2.0, z = 1.25 X 10" 6 in 
normalised coordinates. For the 2D case the z-value is of course 
irrelevant due to the invariance in the z-direction, but we choose 
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Fig. 6. Side views of the field shown in Fig.|4j again at t= 1.05s 
and 50.8s. The collapse of the magnetic field lines in the CMT 
model is obvious. 



it to be small, but non-zero, for the 3D case to avoid creating a 
non-generic orbit. 

T he value for the magn etic moment was also kept the same 
as in IGiuliani et al.1 d2005l) . Because the magnetic fields at the 
starting positions are now different, keeping the magnetic mo- 
ment the same means the i nitial energy of the pa rticles is differ- 
ent to the 6.5keV used by IGiuliani et alj J2005). The values of 
the new magnetic fields at this starting point do not differ signif- 
icantly, so the initial energies are of a similar magnitude to the 
previous work. 

Figure [7] shows the particle orbit for an electron in the 2D 
fields with shear flow. The particle follows the untwisting field- 
lines, and this can be seen clearly in the projections of the trajec- 
tory onto the coordinate planes, which are shown on the sides of 
the box. The orbit looks otherwise similar to th e 2D case without 
shearing as examined by Giuliani e t al.l d2005l) . 

The kinetic energy of the particle as it travels through the 
trap is shown in Fig. [8] As in the 2D case, the energy is gained 
initially mainly due to the effects of the curvature drift, whereas 
in later stages the betatron effect is stronger. The particle starts 
with an energy of 6.5 keV. After 95 seconds the particle energy 
has increased by a facto r of about 6 to 38.0 ke V. This is a simi- 
lar gain to that seen by Giuliani et al. (2005) using the stretched 
field without shear flow to accelerate an electron with initial en- 
ergy of6.5keV to 37.3 keV. 

The particle orbit in the 3D collapsing trap is shown in Fig. [9] 
This shows the effect of the untwisting fieldlines on the particle 
trajectory. A notable difference from the 2D CMT with shear 
flow is the asymmetric projection of the orbit onto the x-z- and 
y-z-planes, wh ereas in the x-y-plan e the orbit looks very similar 
to the orbit in lGiuliani et al.l(l2005l) . 

The energy of the electron in the 3D example is plotted in 
Fig. [10] As the 3D magnetic field decreases faster with height 
than the 2D field, the initial particle energy is lower for the same 
magnetic moment than in the 2D case with shear flow. The initial 
energy is about 3 keV and increases to about 16 keV after 95 
seconds, which corresponds to an increase by a factor 5, whereas 
in the 2D cases we had an increase by a factor of just short of 6. 

We remark that for both examples presented here, we have 
not yet tried to find initial conditions which give rise to higher 
energy gains than rather modest ones that we have found in 
our examples. A more systematic investigation of the 2D CMT 



Grady & Neukirch: 2.5D and 3D Collapsing Magnetic Trap Models 



7 




Fig. 7. Particle orbit in the 2D CMT with shear flow. Projections 
of the trajectory onto the coordinate planes are shown on the 
sides of the box. 
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Fig. 8. Time evolution of the particle energy in the 2D CMT with 
sh ear flow. This evoluti on is very similar to the orbit discussed 
bv lGiuliani etaT] d2005l) . 



of iGiuliani etafl (|2005) shows that much larger energy gains 
are po ssible with increases of a factor 50 or more dGradv et al.l 
2009). We expect similar energy gains to be possible for the 
cases presented here. Another reason for the rather modest in- 
crease in energy is that we have been conservative in our as- 
sumptions about the maximum magnetic field strength on the 
photosphere, which is only about 100 G. A factor 5 to 10 in- 
crease of the photospheric field strength seems reasonable, in 
particular for flaring regions, and this could have a significant 
effect on energy gain. We plan to investigate this in the future. 




Fig. 9. Particle orbit in the 3D CMT model. 
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Fig. 10. Energy gain of the particle in the 3D CMT model. 



5. Summary and Conclusions 

We have developed a fully analytical model for kinematic time- 
dependent, 2D and 3D collapsing magnetic traps. This kinematic 
approach has the advantage that it allows us full control over all 
the features of the model, but has the disadvantage the modelling 
of the plasma system is not self-consistent. 

In the present paper, we have shown how to build kinematic 
CMT models using the magnetic field directly, rather than us- 
ing a flux function or Euler potentials. This is much easier and 
more straightf orward to use, esp ecially in 3D, than the theory 
presented in IGiuliani et al.l d2005l) . 

We have given illustrative examples of collapsing traps with 
transformations that give rise a shear flow in 2D and magnetic 
twist in 3D. We have calculated particle orbits for these new 
CMT models using guiding centre theory. For those orbits, the 
CMT models were found to give similar relative energy gains as 
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with the 2D CMT without shearin g. The particle orbi ts are dif- 
ferent from the 2D CMT model bv lGiuhani et all d2005l) despite 
starting from the same initial position due to the differences in 
field line motion caused by the shear flow in 2D and by the twist- 
ing motion in 3D. The examples shown in this paper have been 
ch osen specifica l ly to b e comparable with the example shown 
in GiulianTetaD(2005). Different CMT models could allow for 
higher energy gains and will be considered in future work. There 
are also many other possible combinations of initial positions, 
initial particle energy and pitch angles, as well as investigating 
proton/ion orbits as well as electron orbits. A sys tematic investi- 
gation for the 2D model of Giuli ani et alj (H005) has shown that 
energy gain factors of or der 50 or higher are possible for that 
model dGradv et al J2 009). A similar investigation is planned for 
the future for 2D with shear flow CMTs and 3D CMTs using the 
theory presented in this paper. 

Acknowledgements. We thank the anonymous referee for useful comments. The 
authors acknowledge financial support by the UK's Science and Technology 
Facilities Council and by the European Commission through the SOLAIRE 
Network (MTRN-CT-2006-035484). 



References 

Aschwanden, M. J. 2002, Space Science Reviews, 101, 1 
Aschwanden, M. J. 2004, ApJ, 608, 554 
Aschwanden, M. J. 2009, Asian J. Phys., 17, 423 
Birn, J., Fletcher, L., Hesse, M., & Neukirch, T. 2009, ApJ, 695, 1151 
Birn, J., Thomsen, M. R, Borovsky, J. E., et al. 1997, J. Geophys. Res., 102, 
2325 

Birn, J., Thomsen, M. R, Borovsky, J. E., et al. 1998, J. Geophys. Res., 103, 
9235 

Birn, J., Thomsen, M. R, & Hesse, M. 2004, Phys. Plasmas, 11, 1825 
Bogachev, S. A. & Somov, B. V. 2001, Astronomy Reports, 45, 157 
Bogachev, S. A. & Somov, B. V. 2005, Astronomy Letters, 31, 537 
Bogachev, S. A. & Somov, B. V. 2009, Astronomy Letters, 35, 57 
Emslie, A. G., Dennis, B. R., Holman, G. D., & Hudson, H. S. 2005, 

J. Geophys. Res., 110, 11103 
Emslie, A. G., Kucharek, H, Dennis, B. R., et al. 2004, J. Geophys. Res., 109, 

10104 

Fletcher, L. & Hudson, H. S. 2008, ApJ, 675, 1645 
Forbes, T. G. & Acton, L. W. 1996, ApJ, 459, 330 
Giuliani, P., Neukirch, T., & Wood, P. 2005, ApJ, 635, 636 
Grady, K. J., Neukirch, T., & Giuliani, P. 2009, in preparation 
Hesse, M, Forbes, T. G., & Birn, J. 2005, ApJ, 631, 1227 
Hesse, M. & Schindler, K. 1988, J. Geophys. Res., 93, 5559 
Hoshino, M., Mukai, T, Terasawa, T, & Shinohara, I. 2001, J. Geophys. Res., 
106, 25979 

Karlicky, M. & Barta, M. 2006, ApJ, 647, 1472 
Karlicky, M. & Kosugi, T. 2004, A&A, 419, 1159 
Kovalev, V. A. & Somov, B. V. 2002, Astronomy Letters, 28, 488 
Kovalev, V. A. & Somov, B. V. 2003a, Astronomy Letters, 29, 111 
Kovalev, V. A. & Somov, B. V. 2003b, Astronomy Letters, 29, 409 
Krucker, S., Battaglia, M., Cargill, P. J., et al. 2008, A&A Rev., 8 
Miller, J. A., Cargill, P. J., Emslie, A. G., et al. 1997, J. Geophys. Res., 102, 
14,631 

Moffatt, H. K. 1978, Magnetic field generation in electrically conducting fluids 

(Cambridge, England, Cambridge University Press, 1978. 353 p.) 
Neukirch, T. 2005, in ESA Special Publication, Vol. 600, The Dynamic Sun: 

Challenges for Theory and Observations 
Neukirch, T, Giuliani, P., & Wood, P. D. 2007, in Reconnection of Magnetic 

Fields, ed. J. Birn & E. Priest (Cambridge University Press), 281-291 
Northrop, T. 1963, The Adiabatic Motion of Charged Particles (New York: 

Interscience Publishers, Inc) 
Piatt, U. & Neukirch, T. 1994, Sol. Phys., 153, 287 

Reeves, K. K, Guild, T. B., Hughes, W. J., et al. 2008a, J. Geophys. Res., 113,0 

Reeves, K. K, Seaton, D. B., & Forbes, T. G. 2008b, ApJ, 675, 868 

Romeou, Z. & Neukirch, T. 1999, in ESA Special Publication, Vol. 448, 

Magnetic Fields and Solar Processes, ed. A. Wilson & et al., 871 — h 
Romeou, Z. & Neukirch, T. 2002, Journal of Atmospheric and Solar-Terrestrial 

Physics, 64, 639 

Somov, B. 1992, Astrophysics and Space Science Library, Vol. 172, Physical 

Processes in Solar Flares (Kluwer Academic Publishers) 
Somov, B. V. & Bogachev, S. A. 2003, Astronomy Letters, 29, 621 



Somov, B. V. & Kosugi, T. 1997, ApJ, 485, 859 

Stern, D. P. 1970, American Journal of Physics, 38, 494 

Stern, D. P. 1987, J. Geophys. Res., 92, 4437 

Veronig, A. M., Karlicky, M., Vrsnak, B., et al. 2006, A&A, 446, 675 



Appendix A: Detailed calculation for the 3D case 
using Euler potentials 

For the following derivation we use a notation which allows us 
to handle as vectors certain groups of scalar quantities or rows or 
columns of tensors. Firstly, the derivatives of the Clebsch vari- 
ables, a and p, with respect to the transformed coordinates are 
required: 



da 
dX 

dp_ 

dX 



da da da 
dX'dY'dZ 
dj3_ dp_ dp 
dX'dY'dZ 



which is basically the usual gradient with respect to X, Y and 
Z. We also need the transformation differentiated with respect to 
the original Eulerian coordinates. 



dX , 


(dX dY dZ\ 


~dH ~ 1 


[dx' dx' dx) 


dX | 


IdX dY dZ\ 


dy ~ 1 


\dy' dy' dy) 


dX | 


(dX dY dZ\ 


~dl ~ ' 


\dz' dz' dz) 



We now consider each component of Eq. (fl5b . starting with the 
x-component 



x dy dz dz dy 



(A.l) 



da 


dX\ 


(dp 


dX\ 


(da 


dX\ 


(dp 


dX 


dX 


dy) 


\dX 


dz) 


[dX 


dz) 


\dX 


dy 



With the coordinate transformation, Eqs.[l8]and ( fl9l l, and using 
the chain rule this becomes 



(A.2) 
(A.3) 
(A.4) 

(A.5) 

(A.6) 
(A.7) 



Applying the well-known vector identity 

(A C)(B D)-(A D)(B ■ O = (A X E) ■ (C X D) 

to Eq. ( IA.2b we arrive at 

•-(£"£)■*»• 

because the initial magnetic field is 

„ _ da dB 
B Q (X) = — x — 

' dx dx 

by construction. Similarly one finds that 
(dX dX\ 



By = 



dz*^-* ™' 



dx dy ) 



